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Abstract 

A calculation of the classical analogue for the quantum wave function and local density of states, in energy 
representation, is presented for simple Hamiltonian systems. Such analogous were proposed by M. V. Berry 
and A. Voros considering the intersection of energy shells of two systems as the only semiclassical object 
which can give support to eigenfunctions. One of them is the system under study and the other one is the 
unperturbed system used to express the wave functions, even in the case that both systems are not close. 
For simple systems and as for scalable ones analytical expressions are obtainable. In the present work we 
offer examples of both. 
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1. Introduction 

Quantum classical correspondence has been an im- 
portant issue since the creation of quantum theory. 
The general idea is expressed in Bohr's principle of 
correspondence. This subject is matter of current re- 
search in a wide class of aspects [Tl. Here we discuss 
the case of eigenfunctions. Several studies and ap- 
proximation to the problem exists in the literature 
y but we shall concentrate in the proposal by M. V. 
Berry and A. Voros Q. The argument is that the 
energy shell of the Hamiltonian systems is the only 
semi classical object who could give support to eigen- 
functions. The practical development of such an idea 
was performed by Benet et al 5] . 

The specific form of such object is discussed in the 
Methodology section. This subject with some correc- 
tions @, 01 , is expected to be the classical equivalent 
of quantum eigenfunctions, and we shall call them 
the Classical Eigenfunctions, CEF. In order to test 
this proposal and due to the nature of this object 
it is expected that a better quantum classical corre- 



spondence occurs in ergodic systems, naturally, the 
easiest way is to consider systems with a chaotically 
classical dynamics. Several systems have been con- 
sidered in order to test the quantum classical corre- 
spondence, and, as a consequence the reliability of the 
CEF (see equation ([1]) below) as the classical equiv- 
alent for the quantum eigenfunctions. The systems 
considered are chaotic quartic oscillators 0, and 
chaotic billiards [Si] , all them with good results. As 
pointed out in [7|, a better fitting of CEF to quan- 
tum eigenfunctions is obtained when both hamiltoni- 
ans, the system under study and that used as a basis, 
present a chaotic classical dynamics. This is a result 
that CEF depends on the ergodic properties of the 
systems. But what about integrable systems, how 
good could be the quantum classical correspondence 
in these systems?. Another important point is that in 
all previous works the CEFs were calculated numer- 
ically testing carefully the reliability, however no so 
much was clarified about the properties of them. An 
example of this is that in all the reported cases CEFs 
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shown strong peaks, but nothing is discussed about 
their nature. Here, we start a small contribution on 
the subject. 

2. Methodology 

Looking for a practical formulation of Berry's and 
Voros' idea, L. Benet etal 5] proposed that the quan- 
tity, 

g{£, £o)^aJ 5{£ - H)S{£o - Ho)dpdq, (1) 

could be the classical counterpart for the eigenfunc- 
tions of the hamiltonian H at energy £ in the basis of 
the hamiltonian Hq with energies £q. A is a normal- 
ization constant which contains information about 
the quantum scale occupied by the state in phase 
space and 6 is the Dirac's delta function. This quan- 
tity does not contain information about the quan- 
tum phase, neither is the solution of an eigenfunction 
equation. However, it corresponds to the intensity of 
a quantum eigenfunction j^*) in the |^o) basis; for 

H\^)=£\^), (2) 

and 

HolM =£o\<f>o), (3) 

respectively. Evidence of this quantum classical cor- 
respondence for hamiltonian systems with a chaotic 
classical counterpart exists The function g 

represents the intersection of both energy shells, the 
H and Hq ones, and it is obtained in the same way 
as the average density of states though the Weyl for- 
mula (see, for instance, Q). It is important to recall 
that distributions like (P) or the Weyl approximation 
requires certain ergodicity in the system. 

A way to understand the physical meaning of ([Ij 
go as follow: Consider a system H composed by the 
sum of Hq and a perturbation XW, i.e., H — H+XW. 
The perturbation is not necessarily small. At a fixed 
energy £ of the system H, the energy is distributed 
between both parts, the unperturbed and the pertur- 
bation, hence the energy £q is time dependent and 
its distribution defines ([T]). In other words, if we con- 
sider the classical solution of the test system and we 
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Figure 1: (a) Time evolution of unperturbed energy at fixed 
system energy, (b) Histogram of unperturbed energy distri- 
bution for the case sliown in (a). In log-linear scale. At the 
mixing time the distribution becomes JlJ. Units are arbitrary. 

evaluate these solutions in phase space in the non- 
perturbed hamiltonian, the corresponding energy is 
a time dependent variable. In Fig. la we show this 
process. As the time pass, a limit distribution ap- 
pears as is shown in Fig. lb. The distribution is the 
histogram of the projection into the unperturbed en- 
ergy for a time sufficiently large. The time needed to 
reach such a distribution is called the mixing time. 
Meanwhile this picture is clear, a practical calcula- 
tion of the distribution depends strongly on the ini- 
tial conditions. Here and in 0, 0, 0] the function in 
(HI) is calculated using Montecarlo integration. This 
method allows an efficient calculation. 

2.1. A small but important detail 

Notwithstanding that ([1]) shows a good quantum 
classical correspondence in the cases presented in 
@, Sj; expression ([T|) must be rectified in order to 
have a proper normalization. As we shall see, nor- 
malization depends on both energies £, and £o via 
their corresponding density of states. This detail be- 
comes much more important as we consider a larger 
number of particles and the densities become larger 
as well. The definition proposed [1, 0] is 

gc{£,£o) = {c/p{£)pq{£q)) 

X /^(£-ff)<5(£o-i?o)dpdq, 



where p{£) and po{£o) are the density of states for 
the perturbed and unperturbed systems respectively. 
They are calculated according to the Weyl approxi- 
mation 

pix) = (1/c) J 5{x - i/)dpdq, (5) 

for any of the hamiltonians. The constant c is de- 
fined as l/{2TihY. Here, d is the number of degrees 
of freedom of our systems. Even when calculations 
could be performed in three dimensional configura- 
tion space for any number of particles we shall con- 
centrate in one dimensional systems and hence d will 
be the number of particles considered. 

Now we shall explain how to arrive to correction 
(HI) from ([T]). The arguments in order to arrive to 
expression (jj]) are as follow @: 

Integrating ^ respect to £o we obtain 

J g{£,£o)d£o^ p{£){27Thf (6) 

where we used the properties of the delta function. 
An integration on £ gives as a result po{£o{27rh)'^). 
The generalization to the continuum of the normal- 
ization rule of a quantum state on a discrete or a 
continuum spectrum requires an integration on the 
energy element po{£o)d£o. Hence, only ^ offers the 
correct normalized function that fulfills all these re- 
quirements. 

For practical porpoises we require a simplified ver- 
sion of ([4]). For simplicity we shall consider systems 
of one or several particles in one dimension, hence 
the parameters d shall refer to the number of par- 
ticles. In this section, as well, we shall refer to ([Ij 
instead of ([4]) in order to avoid to carry on the densi- 
ties. In absence of magnetic filed, the integration of 
the momenta variables can be performed easily with 
the help of the property 

Six - mx -g) = 6{x - (/ + g)l2)5{J - <?), (7) 

for any scalar functions / and g. This allows to 
rewrite the integrand of ([iJ as 

5{p^/2-{£-V+£o-Vo)/2)5{£-V-{£a-V(,)), (8) 



for the potentials associated to each hamiltonian. Us- 
ing spherical coordinates for the momenta in the ar- 
gument of ([5]), the integral defining (|H) can be rear- 
ranged as 

g,{£,£o) = B^dq_5{£-V -{£o-Vo)) 

d-2 \ 

x(£-V + £Q-Va))—Q{{£-V + £Q-V^)). 

Here is the step function and appears in order 
to keep the positivity in the square root argument. 
B stands for 

B = 27:"/^/ {{2^hfV{d/2)p(£)po{£a)) (10) 

Where F stands for the Gamma function. This con- 
stant contains the result from the angular part in the 
momentum integral. We shall use extensively ^ for 
the calculations presented in this work. 

Property ([7]) can be obtained straightforwardly us- 
ing the sequence of gaussians with limit a Dirac delta 
function. Assuming that the sequence if n gaussians 
is taken at the same time for each case, we can ar- 
range properly the argument in both gaussians and 
take the limit to infinity in order to obtain the deltas 
again(see appendix in [6|). 

2.2. The systems 

In order to understand the use and the properties 
of ^ we shall consider several simple systems and 
one not so simple in order to contrast our results. 
The selected systems are a general one dimensional 
Hamiltonian in their own basis, the free particle in its 
own basis, a quartic oscillator in the harmonic basis 
and, at the end, a chaotic several particles Hamilto- 
nian in the basis of the quartic single particle basis. 

3. Results 

3.1. Simple special cases in one dimension 

Here we apply the considerations discussed in pre- 
vious section. We discuss first some very simple ex- 
amples in one dimension for one particle. In subsec- 
tion B we shall discuss the cases with more than one 
particle. 
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3.1.1. Equal Hamiltonians in d dimensions 
Eigenfunctions of quantum systems written in their 

own basis must be Kronecker's or Dirac's deltas. This 
is the case for the classical counterpart as well. To 
show this we use ([T]) and we consider as an indepen- 
dent variable the whole Hamiltonian H, hence we 
have 

g{E,Eo) = /dqdq,5(f -F5(fo-i?o)), , . 

= 6{E^Eo){2TTh)ip{E). ^ ' 

We used property ^ in order to get the last line. 
The CEF could be written as 

gc{£,£o) = 5{£-£a)p^{£o). (12) 

Notice that, as expected, these results do not depend 
on the potentials used or the system's dimensionality. 

3.1.2. Free particle in its own basis 

This a particularly simple case, the free particle 
in one dimension with H = Hq = p'^/2. Using the 
previous dimension with result (|12p we only require 
to evaluate the density of states pq{£q), giving 

po{£o) = il/27rh)JSipy2~£o)dpdq 

= {l/27Th){L/V2£). ^ ' 

Where we performed the required integral in g on a 
box of size L and we take the limit to infinity at the 
end. The final result is 

£o) = {27rhV2£)Si£ - £o). (14) 

The dependence on L disappears when the normal- 
izations is performed. The same result is obtained by 
direct integration. 

3.1.3. Quartic oscillator in the harmonic basis. 
Now we consider the perturbed Hamiltonian in the 

H ^ /2 + Pq^ + aq'^ basis of the harmonic hamil- 
tonian, i.e., H — /2 + Pq^. Using ([SJ with d 1 to 
this case we have 



5c (X y dq5{f{q))/^£ - (/3q2 + aq^) + (£„ _ pq2)^ 

(15) 



with 

f{q)=£-{pq^+aq^)-{£^~pq^), (16) 

and we dropped the unit step function for sake of 
clarity. Eq. (jl5p can be re-written as 



ffc (X y dq5{£ ~£o- aq^)/^£ + £o-i2f3q^+aq^). 

(17) 

We require the roots in the argument of the delta 
function, that is 

q = {/{£ - £o)/a. (18) 
Using the property of delta functions 

^(/(a;))-E'^(^-^^)/l/'(^')l' (19) 

i 

on (IT7| . we obtain 

cx {{£ - £a)/a)-^'^ (fo - N{£ - fo)/a) . 

(20) 

The density of states is constant for both Hamil- 
tonians, the unperturbed one is proportional to the 
mean level density, which, in this case, is constant. 
(The spectrum of the one dimensional harmonic os- 
cillator is equidistant). The density of states for the 
perturbed case can be calculated explicitly with the 
help of symbolic calculus software, but at the end 
it have a is a number. From here it is clear that 
gc((£,£Q) have a divergence when both energies are 
equal. Notice, as well, that another classically forbid- 
den region appears when the argument the asymme- 
try in the root argument £o — P^J {£ — £o)/a becomes 
negative, that is, it corresponds to an asymmetrical 
function. In Fig. 2 we show a graph of this case. 
In broken red line appears the classical analogue (I^IH) 
for the energy corresponding to the state number 756 
of the perturbed system. In black points the corre- 
sponding quantum result. An average on windows of 
size 30 was performed in order to smooth the quan- 
tum oscillations, the result is shown in blue crosses. 
The agreement is notable. A shift (not shown) of the 
classical eigenfunction is near of the quantum case 
envelope. 
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Figure 2: (Color on-line) Eigcnfunction intensities in semi log 
scale for the state 756 of the single particle quartic oscillator in 
the harmonic basis ( black pints), classical equivalent for the 
same case (broken red line). We add an averaged version of 
the quantum eigenfunction. 




[5 ' leO ISO 200 220 2+0 

Unperturbed energy 




Unperturbed energy 



Figure 3: (a) Average intensities for the state number 500 
over 101 eigenfunctions around it (fluctuating curve) and their 
corresponding classical analogous (smooth curve), (b) Same 
as before but in semi log scale. 



Notice that the non-oscillating tails in the quan- 
tum case, shown in Fig 2, correspond to the tunnel 
effect to the classically forbidden zone. The quantum 
case was calculated using the diagonalization of the 
Hamiltonian matrix in the unperturbed basis using 
4000 states. 

3.2. Not so simple cases 

Analogues like ([9]) are builded up for complex sys- 
tems where a certain amount of ergodicity and per- 
haps mixing exists. For such cases quantum classical 
correspondence has been good 0, H, 0, II], here we 
shown another example. The Hamiltonian consid- 
ered is equal to those appearing in 0] 

(21) 

with the following values for the parameters a = 10.0, 
/? = —2.1 and 7 = —2.2, and the number of particles 
is d = 4. Here, T stands for the kinetic energy. The 
numerical calculation was performed diagonalizing a 
10000 X 10000 matrix with the described method in 
[lo| and a proper basis cutoff. The chaoticity was 
tested calculating the maximum Lyapunov exponent 
for a large variety of initial conditions. In Fig. 3 we 
plot an average of 100 eigenfunctions around the state 
number 500. The quantum classical correspondence 



is good but a shoulder of quantum nature appears. 
This shoulder causes that the fit at the tail fails as a 
result of the normalization rule. Notice that the clas- 
sical analogue present a strong peak at the resonance 
energy meanwhile the quantum counterpart present a 
large intensity at the surrounding but not precisely at 
that energy. In fact, the average of quantum eigen- 
functions was performed shifting the eigenfunctions 
in order that the maximum intensity coincide in each 
case. This is possible due that we are considering 
eigenstates sufHciently high in the energy spectrum. 

4. Discussion 

Even when the reliability of the numerical calcula- 
tions is achieved, several characteristics of CEF (|9]) 
remains obscure, that is the motivation of this simple 
calculations. As explained in previous sections, the 
analytical results depends in depth on the nature of 
the potentials used. However some general features 
could be discussed. First, the square root and the 
positivity of the arguments is of general nature for 
conservative systems, hence such characteristics will 
remain even when we use several particle Hamiltoni- 
ans in three dimensions. These characteristics could 
be important since eigenfunctions equivalents could 
give a faster estimation for some wave functions of 
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interest in atomic or molecular physics. Addition- 
ally, the zeros of the arguments in Q or (|9l) define 
the border to the classically forbidden zone and the 
turning points. Second, one of main problems in or- 
der to obtain completely the integrals involved in ^ 
is that we require the knowledge of the zeros of a 
complicated multidimensional function, and express 
them in a generalization of (|18l) . 

Another interesting point to discuss is about the 
CEF behavior £q near £. From the numerical cal- 
culations this behavior at point seems like a diver- 
gence and not a smooth peak. The examples pre- 
sented here have a divergence at this point. Hence 
for all these cases a divergence peakshould appears 
in the case that the unperturbed and the perturbed 
energies coincide. Others divergences associated with 
classical turning points must appear in the evaluation 
of dH) as appears in the case presented in 13.1.31 How 
CEF looks like for a larger number of particles in this 
case is not obvious since degenerancies appear in the 
several particles spectrum, which invites to look for 
another example where the generalization to several 
particles and simplicity of analytical calculations can 
be applied. 

For a large number of particles, the limit expres- 
sion of dl]) tend to be a gaussian like form [6]. This 
shape is the expected one for general eigenfunctions 
like those generated by gaussian ensembles of matri- 
How eigenfunctions change this behavior is 
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ces 

not clear, but even for the numerical calculations pre- 
sented here and those in progress for a larger number 
of particles the peaks appear as a divergence. 

Up to now we deal with eigenfunctions, that is, 
we consider CEF Q and ([9]) as a function of unper- 
turbed energy at fixed £ . For the local density of 
states the oposit is considered, (|4]) and dH) become 
functions of £ instead of £q. This case is of relevance 
for experimental comparison. Due to the symmetry 
in (|4]) we can use, mutatis mutandis, the same expres- 
sions as before, including In Fig. 4 we shown 
the LDOS calculation for the same case presented 
in section 13.21 In smooth line the classical counter- 
part obtained from (U) and in strongly oscillating line 
the corresponding quantum case averaged in a similar 
way as before. Notice that the quantum classical cor- 
respondence is much better fitted, and the shoulder 
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Figure 4: Local Density of States for the anharmonic oscillator 
presented in the text. Classical analogue in smooth curve. 



is not present. 

5. Conclusions 

Classical analogues for eigenfunctions. Local Den- 
sity of States and occupation number statistics have 
been proposed in 0, [a, 0, S] with a good quantum 
classical correspondence for few particle anharmonic 
hamiltonian and for a chaotic billiard. Evaluation of 
such analogues has been performed numerically and 
some characteristics remain obscure. Here we started 
a program in order to understand general properties 
of those quantities, like divergences, ergodicity and 
tunnel effect. For that matter we evaluate analyti- 
cally the expression for the classical equivalents for 
very simple cases. The cases were selected for sim- 
plicity and in order to identify the causes of diver- 
gence without the use of a numerical calculation. The 
cases presented here shown that the numerical eigen- 
functions calculated have a divergence when the en- 
ergy coincides with the perturbed one. An interesting 
point is the limit where a large number of particles 
is considered. In such a case eigenfunctions shown a 
rounded peak at the resonance energy and a gaussian 
form, unfortunately the harmonic oscillators consid- 
ered in the present work have a lot of degenerancies 
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and them avoid to understand the hmit, even when 
in the case presented here have a notable good quan- 
tum classical correspondence for one single particle 
and without the requirement of ergodicity. 

As a final remark we should stress that classi- 
cal eigenfunction analogue (|4]) is well situated for a 
fast calculation of average eigenfunctions in quantum 
complex systems. 
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